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(54) Method and system for multi-processor FFT/IFFT with minimum inter-processor data 
communication 



(57) The present invention provides a scalable 
method for implementing FFT/IFFT computations in 
multiprocessor architectures that provides improved 
throughput by eliminating the need for inter-processor 
communication after the computation of the first "log 2 P" 
stages for an implementation using "P" processing ele- 
ments, comprising computing each butterfly of the first 
"log 2 P" stages on either a single processor or each of 
the "P" processors simultaneously and distributing the 
computation of the butterflies in all the subsequent stag- 



es among the "P" processors such that each chain of 
cascaded butterflies consisting of those butterflies that 
have inputs and outputs connected together, are proc- 
essed by the same processor. 

The invention also provides a system for obtaining 
scalable implementation of FFT/IFFT computations in 
multiprocessor architectures that provides improved 
throughput by eliminating the need for inter-processor 
communication after the computation of the first "log 2 P" 
stages for an implementation using "P" processing ele- 
ments. 
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Description 

[0001] The present invention relates to the field of dig- 
ital signal processing. More particularly the invention re- 
lates to a device and method for providing an FFT/IFFT 
implementation providing minimum inter-processor 
communication overhead and less silicon area in a mul- 
tiprocessor system. 

[0002] The class of Fourier transforms that refer to 
signals that are discrete and periodic in nature are 
known as Discrete Fourier Transforms (DFT). The dis- 
crete Fourier transform (DFT) plays a key role in digital 
signal processing in areas such as spectral analysis, fre- 
quency domain filtering and polyphase transformations. 
[0003] The DFT of a sequence of length N can be de- 
composed into successively smaller DFTs. The manner 
in which this principle is implemented falls into two class- 
es. The first class is called a "decimation in time" ap- 
proach and the second is called a "decimation in fre- 
quency" method. The first derives its name from the fact 
that in the process of arranging the computation into 
smaller transformations the sequence "x(n) n (the index 
'n 1 is often associated with time) is decomposed into suc- 
cessively smaller subsequences. In the second general 
class the sequence of DFT coefficients "x(k)" is decom- 
posed into smaller subsequences (k denoting frequen- 
cy). The present invention employs "decimation in time". 
[0004] Since the amount of storing and processing of 
data in numerical computation algorithms is proportional 
to the number of arithmetic operations, it is generally 
accepted that a meaningful measure of complexity, or 
of the time required to implement a computational algo- 
rithm, is the number of multiplications and additions re- 
quired. The direct computation of the DFT requires 
"4N 2 " real multiplications and "N(4N-2) N real additions. 
Since the amount of computation and thus the compu- 
tation time is approximately proportional to "N 2 " it is ev- 
ident that the number of arithmetic operations required 
to compute the DFT by the direct method becomes very 
large for large values of "N". For this reason, computa- 
tional procedures that reduce the number of multiplica- 
tions and additions are of considerable interest. The 
Fast Fourier Transform (FFT) is an efficient algorithm 
for computing the DFT 

[0005] The conventional method of implementing an 
FFT or Inverse Fourier Transform (I FFT) uses a radix-2 
/ radix-4 / mixed-radix approach with either "decimation 
in time (DIT)" or a "decimation in frequency (DIF)" ap- 
proach. 

[0006] The basic computational block is called a "but- 
terfly" - a name derived from the appearance of flow of 
the computations involved in it. Fig. 1 shows a typical 
radix-2 butterfly computation. 1.1 represents the 2 in- 
puts (referred to as the "odd" and 'even" inputs) of the 
butterfly and 1 .2 refers to the 2 outputs. One of the inputs 
(in this case the odd input) is multiplied by a complex 
quantity called the twiddle factor (W N k ). The general 
equations describing the relationship between inputs 



and outputs are as follows: 

X[k] = x[n] + x[n+N/2]W N k 

5 

X[k+N/2] = x[n] - x[n+N/2]W N k 

[0007] An FFT butterfly calculation is implemented by 
10 a z-point data operation wherein "z" is referred to as the 
"radix". An "N" point FFT employs "N/z" butterfly units 
per stage (block) for "log z N" stages. The result of one 
butterfly stage is applied as an input to one or more sub- 
sequent butterfly stages. 
15 [0008] Computational complexity for an N-point FFT 
calculation using the radix-2 approach = 0(N/2 * log 2 N) 
where "N" is the length of the transform. There are ex- 
actly "N/2 * log 2 N" butterfly computations, each com- 
prising 3 complex loads, 1 complex multiply, 2 complex 
20 adds and 2 complex stores. A full radix-4 implementa- 
tion on the other hand requires several complex load/ 
store operations. Since only 1 store operation and 1 load 
operation are allowed per bundle of a typical VLIW proc- 
essor that is normally used for such implementations, 
25 cycles are wasted in doing only load/store operations, 
thus reducing ILP (Instruction Level parallelism). The 
conventional nested loop approach requires a high loop- 
ing overhead on the processor. It also makes application 
of standard optimization methods difficult. Due to the na- 
30 ture of the data dependencies of the conventional FFT/ 
IFFT implementations, multi cluster processor configu- 
rations do not provide much benefit in terms of compu- 
tational cycles. While the complex calculations are re- 
duced in number, the time taken on a normal processor 
35 can still be quite large. It is therefore necessary in many 
applications requiring high-speed or real-time response 
to resortto multiprocessing in order to reduce the overall 
computation time. For efficient operation, it is desirable 
to have the computation as linearly scalable as possible 
40 - in other words the computation time reducing in inverse 
proportion to the number of processors in the multiproc- 
essing solution. Current multiprocessing implementa- 
tions of FFT/IFFT however, do not provide such a linear 
scalability. 

45 [0009] US patent 6,366,936 describes a multiproces- 
sor approach for efficient FFT. The approach defined is 
a pipelined process wherein each processor is depend- 
ent on the output of the preceding processor in order to 
perform its share of work. The increase in throughput 
50 does not scale proportionately to the number of proces- 
sors employed in the operation. 
[0010] US patent 5,293,330 describes a pipelined 
processor for mixed size FFT. Here too. the approach 
does not provide proportional scalability in throughput, 
55 as it is pipelined. 

[0011] A scheme for parallel FFT/IFFT as described 
in "Parallel 1-D FFT Implementation with TMS320C4x 
DSPs" by the semiconductor group Texas Instruments, 
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uses butterflies that are distributed between two proc- 
essors. In this implementation, inter-processor commu- 
nication is required because subsequent computations 
on one processor depend on intermediate results from 
other processors. Every processor computes a butterfly 
operation on each of the butterfly pairs allocated to it 
and then sends half of its computed result to the proc- 
essor that needs it for the next computation step and 
then waits for the information of the same length from 
another node to arrive before continuing computation. 
This interdependence of processors for a single butterfly 
computation does not support proportionate increase in 
output with increase in the number of processors. 
[0012] European patent application no. 03027181.1 
filed on 27.11 .03 describes a linearly scalable FFT/IFFT 
system. The system incorporates a shared memory 
wherein each processor accesses correct data samples 
from the shared memory. Distribution is such that no in- 
ter-processor communication is required during the but- 
terfly computation. However there is a requirement of 
inter-processor communication between stages. 
[0013] Though a shared memory system is easier it 
is not very economical. This is because this approach 
requires multiport memories that are very expensive. 
Therefore a distributed memory system is more eco- 
nomical. The distributed memory architecture requires 
a media to communicate data among the processors. 
Hence it is desirable that the data communication 
among the processors in distributed memory architec- 
ture is minimum. Since the input data is distributed in 
equal size segments to each processor and each proc- 
essor performs computations only on the data in its local 
memory, the memory requirement for individual proces- 
sor reduces resulting in a lower requirement for silicon 
area and cost. 

[0014] The aim of the present invention is to over- 
come the above drawbacks and provide a device and 
method for implementing FFT/IFFT with minimum com- 
munication overhead among processors in a multiproc- 
essor system using distributed memory. 
[0015] According to the invention, a scalable method 
for implementing FFT/IFFT computations in multiproc- 
essor architectures is provided, according to claim 1 . 
[0016] In practice, the present method provides im- 
proved throughput by eliminating the need for inter- 
processor communication after the computation of the 
first "log 2 P" stages for an implementation using "P" 
processing elements. Specifically, the method compris- 
es the steps of: 

computing each butterfly of the first "log 2 P" stages 
on either a single processor or each of the "P" proc- 
essors simultaneously, 

distributing the computation of the butterflies in all 
the subsequent stages among the "P" processors 
such that each chain of cascaded butterflies con- 
sisting of those butterflies that have inputs and out- 
puts connected together, are processed by the 



same processor. 

[001 7] The distributing of the computation of the but- 
terflies subsequent to the first "log 2 P" butterflies is 

5 achieved by assigning operand addresses of each set 
of butterfly operands to each processor in such a man- 
ner that the butterfly is processed by the same proces- 
sor that computed the connected butterfly of the previ- 
ous stage in the same chain of butterflies. 

10 [0018] The desired assignment of operand addresses 
is achieved by deriving the address of the first operand 
in the operand pair corresponding to the "i th " stage of 
the computation from the address of the corresponding 
operand in the previous stage by inserting a "0" in the " 

15 (i+1 )th» bit position of the address, while the address of 
the second operand is derived by inserting a "1" in the 
"(i+1) th " bit position of the operand address. 
[001 9] The above method further includes computing 
of twiddle factors for the butterfly computations at each 

20 processor by initializing a counter and then increment- 
ing it by a value corresponding to the number of proc- 
essors "P" and appending the result with a specified 
number of "0"s. 

[0020] The present invention also provides a system 
25 for obtaining scalable implementation of FFT/IFFT com- 
putations in multiprocessor architectures, according to 
claim 5. In detail the system comprises: 

means for computing each butterfly of the first 
30 M log 2 P" stages on either a single processor or each 
of the "P" processors simultaneously, 
addressing means for distributing the computation 
of the butterflies in all the subsequent stages among 
the "P" processors such that each chain of cascad- 
es ed butterflies consisting of those butterflies that 
have inputs and outputs connected together, are 
processed by the same processor. 

[0021] The addressing means comprises addresses 
40 generation means for generating the operand address- 
es of the butterflies subsequent to the first "\og 2 P" but- 
terflies in such a manner that the butterfly is processed 
by the same processor that computed the connected 
butterfly of the previous stage in the same chain of but- 
45 terflies. 

[0022] The address generation means is a computing 
mechanism for generating the address of the first oper- 
and in the operand pair corresponding to the "i th " stage 
of the computation from the address of the correspond- 
50 ing operand in the previous stage by inserting a "0" in 
the "(i+1) th " bit position of the address, and generating 
the address of the second operand by inserting a "1 " in 
the "(i+1) th " bit position of the operand address. 
[0023] The above system further includes a comput- 
es jng mechanism for address generation of twiddle factors 
for each butterfly on the corresponding processor. 
[0024] The present invention will now be explained 
with reference to the accompanying drawings, which are 
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given only by way of illustration and are not limiting for 
the present invention. 

[0025] FIG. 1 shows the basic structure of the signal 
flow in a radix-2 butterfly computation for a discrete Fou- 
rier transform. 

[0026] FIG. 2 shows a 2-processor implementation of 
butterflies for a 16-point FFT, in accordance with the 
present invention. 

[0027] FIG. 3 shows a 4-processor implementation of 
butterflies for a 16-point FFT, in accordance with the 
present invention. 

[0028] Fig 2 shows the implementation for a 1 6 point 
FFT in a 2-processor architecture using the present in- 
vention. Solid lines are computed in one processor, and 
dashed lines in the other. The computational blocks are 
represented by '0'. The left side of each computational 
block is its input (the time domain samples) while the 
right side is its output (transformed samples). The 
present invention uses a mixed radix approach with dec- 
imation in time. The first two stages of the radix-2 FFT/ 
I FFT are computed as a single radix-4 stage. As these 
stages contain only load/stores and add/subtract oper- 
ations there is no need for multiplication. This leads to 
reduced time for FFT/I FFT computation as compared to 
that with full radix-2 implementation. The next stages 
have been implemented as radix-2. The three main 
nested loops of conventional implementations have 
been fused into a single loop which iterates "N/2*(log 2 
N-2))/(n umber of processor)" times. Each processor is 
used to compute one butterfly in one loop iteration. 
Since there is no data dependency between different 
butterflies in this algorithm, both during and between 
stages, the computational load can be equally divided 
among the different processors, leading to a nearly lin- 
ear scalable system. There is no data dependency be- 
tween stages and therefore each processor is able to 
perform the butterfly computations on the data assigned 
to it without communicating with the other processors. 
[0029] The mechanism for assigning the butterflies in 
this manner consists of generating the addresses of in- 
puts such that each processor computes a complete se- 
quence of cascaded butterflies. An N-bit counter, where 
"N" is the number of stages is used to derive the ad- 
dresses used for variables corresponding to each but- 
terfly stage in the computation. Two inputs are generat- 
ed for the two operands of the butterfly. Introducing '0' 
in a specified position of counter generates address for 
input 1. Introducing 1 in a specified position of counter 
generates address for input 2. For address generation 
of twiddle factors a separate counter with a specified 
number of bits is initialized on each processor. The 
counter value is then appended with a specified number 
of zeroes. The counter is incremented by a value corre- 
sponding to the number of processors and appended 
with a specified number of zeroes to get the twiddle fac- 
tor address of the next butterfly stage. 
[0030] Distribution of data samples to the other proc- 
essor can be after stage 1 at the earliest. But to save on 



unnecessary multiplications it can also be done after 
stage 2. No inter-processor communication is required 
once data is distributed. As a result, the red line outputs 
need to collected by one of the processor at the end of 
5 the computation. 

[0031] Fig 3 shows a 4-processor implementation for 
the 1 6-point FFT using this invention. Different linetypes 
represent computation in each of the 4 processors. 
[0032] In the present implementation of the invention 
each processor comprises one or more ALUs (Arithme- 
tic Logic unit), multiplier units, data cache ; and load/ 
store units. Each processor has an individual memory 
and the distribution of butterflies is such that there is no 
inter-processor communication required after the distri- 
bution of data. The distribution of data takes place after 
"log 2 P" stages where "P" is the number of processors. 
Inter-processor communication takes place only before 
and after all the computations have been completed. 
The amount of data communication overhead does not 
increase with an increase in the number of processors. 
[0033] The invention, though described for distributed 
memory can be applied to shared memory systems al- 
so. European patent application no. 03027181.1 cited 
above provides a linearly scalable system for shared 
memory systems only. 
[0034] It will be apparent to those with ordinary skill in 
the art that the foregoing is merely illustrative intended 
to be exhaustive or limiting, having been presented by 
way of example only and that various modifications can 
be made within the scope of the above invention. 
[0035] Accordingly, this invention is not to be consid- 
ered limited to the specific examples chosen for purpos- 
es of disclosure, but rather to cover all changes and 
modifications, which do not constitute departures from 
the permissible scope of the present invention. The in- 
vention is therefore not limited by the description con- 
tained herein or by the drawings, but only by the claims. 



1. A scalable methodfor implementing FFT/I FFT com- 
putations in multiprocessor architectures that pro- 
vides improved throughput by eliminating the need 
45 for inter-processor communication after the compu- 
tation of the first "log 2 P" stages for an implementa- 
tion using "P" processing elements, comprising the 
steps of : 

50 - computing each butterfly of the first "log 2 P" 
stages on either a single processor or each of 
the "P" processors simultaneously, 
distributing the computation of the butterflies in 
all the subsequent stages among the "P" proce 
55 ssors such that each chain of cascaded butter- 

flies consisting of those butterflies that have in- 
puts and outputs connected together, are proc- 
essed by the same processor. 
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2. A method as claimed in claim 1 wherein the distrib- 
uting of the computation of the butterflies subse- 
quent to the first "log 2 P" butterflies is achieved by 
assigning operand addresses of each set of butter- 
fly operands to each processor in such a manner 
that the butterfly is processed by the same proces- 
sor that computed the connected butterfly of the 
previous stage in the same chain of butterflies. 

3. A method as claimed in claim 2 wherein the desired 
assignment of operand addresses is achieved by 
deriving the address of the first operand in the op- 
erand pair corresponding to the "jth" stage of the 
computation from the address of the corresponding 
operand in the previous stage by inserting a "0" in 
the "(i+1) th " bit position of the address, while the 
address of the second operand is derived by insert- 
ing a "1" in the "(i+1) th " bit position of the operand 
address. 

4. A method as claimed in claim 1 further including the 
computing of twiddle factors for the butterfly com- 
putations at each processor by initializing a counter 
and then incrementing it by a value corresponding 
to the number of processors "P" and appending the 
result with a specified number of "0"s. 

5. A system for obtaining scalable implement ation of 
FFT/IFFT computations in multiprocessor architec- 
tures that provides improved throughput by elimi- 
nating the need for inter-processor communication 
after the computation of the first "log 2 P" stages for 
an implementation using "P" processing elements, 
comprising : 

a means for computing each butterfly ofthefirst 
"log 2 P" stages on either a single processor or 
each of the "P" processors simultaneously, 
an addressing means for distributing the com- 
putation of the butterflies in all the subsequent 
stages among the "P" processors such that 
each chain of cascaded butterflies consisting 
of those butterflies that have inputs and outputs 
connected together, are processed by the 
same processor. 

6. A system as claimed in claim 5 wherein the ad- 
dressing means comprises addresses generation 
means for deriving the operand addresses of the 
butterflies subsequent to the first "log 2 P" butterflies 
in such a manner that the butterfly is processed by 
the same processor that computed the connected 
butterfly of the previous stage in the same chain of 
butterflies. 

7. A system as claimed in claim 6 wherein the address 
generation means is a computing mechanism for 
deriving the address of the first operand in the op- 



erand pair corresponding to the "i th " stage of the 
computation from the address of the corresponding 
operand in the previous stage by inserting a "0" in 
the ) th " bit position of the address, and deriving 
5 the address of the second operand by inserting a 
"1" in the "(i+1 ) th " bit position of the operand ad- 
dress. 

8. A system as claimed in claim 5 further including a 
10 computing mechanism for address generation of 
twiddle factors for each butterfly on the correspond- 
ing processor. 
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